Entanglement dynamics via semiclassical propagators in systems of two spins 
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We analyze the dynamical generation of entanglement in systems of two interacting spins initially 
prepared in a product of spin coherent states. For arbitrary time-independent Hamiltonians, we de- 
rive a semiclassical expression for the purity of the reduced density matrix as function of time. The 
final formula, subsidiary to the linear entropy, shows that the short-time dynamics of entanglement 
depends exclusively on the stability of trajectories governed by the underlying classical Hamilto- 
nian. Also, this semiclassical measure is shown to reproduce the general properties of its quantum 
counterpart and give the expected result in the large spin limit. The accuracy of the semiclassical 
formula is further illustrated in a problem of phase exchange for two particles of spin j. 
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I. INTRODUCTION 

When two initially separated quantum systems are led 
to interact with each other they lose their individuality. 
This means that it is no longer possible to express the 
state of one of the systems separately from the other, 
i.e., they have got entangled. The relevance of these 
quantum correlations, which was recognized already in 
the early days of the Quantum Theory, nowadays dis- 
penses with further highlights. Entanglement has defi- 
nitely achieved a prominent place within the quantum 
phenomenology P, 

In a less consensual scenario, foundational questions 
have been posed which try to decipher if and how en- 
tanglement manifests in the classical limit. Surprisingly, 
even though it is hard to conceive a classical image of 
entanglement at a first sight, there exists a number of 
works reporting on the persistence of entanglement in 
semiclassical regime. Although these works agree on this 
essential point, their approaches are clearly different in 
methodology, interpretation, and even on the very notion 
of semiclassical limit. 

In a seminal work Q , Furuya and co-authors have nu- 
merically shown that in the short-time regime entangle- 
ment behaves in accordance with the underlying classical 
dynamics, with accentuated differences between chaotic 
and regular initial conditions. A key ingredient in this 
approach is the use of coherent states, which are used 
as initial states for the dynamics as well as to furnish, 
through a well-defined prescription, the corresponding 
classical structure. The approaches of Refs. 0, Q follow 
the same essence, though the last one focuses on systems 
of two spins. References on the other hand, pro- 

pose to link entanglement with entropic measures defined 
within classical-statistical theories. Still, some authors 
have investigated the semiclassical limit of entanglement 
(and of decoherence) by applying time-dependent pertur- 
bation theory and diagonal approximations |lll4l3l | . 

The present work lies in the context delineated by 
Refs. |14| - [la |. Basically, these papers employ semiclas- 
sical propagators to analyze the entanglement dynam- 



ics of bipartite quantum systems. In Jacquod's ap- 
proach [3, the approximation is performed using 
momentum and space representations simultaneously, 
while in our previous article [l6| we adopt the coherent- 
state representation. Although both calculations yield 
the same basic results, ours has the advantage of having 
been naturally structured to accommodate spin degrees 
of freedom. The aim of this contribution is to carry on 
this program, providing, for the first time, a semiclas- 
sical expression for entanglement dynamics of two-spin 
systems. 

This paper is organized as follows. We start Sec. [Tl] 
by reviewing the main elements of the formal structure 
associated to the semiclassical spin-coherent-state prop- 
agator. We then introduce the time-reversal propagator 
and unify the formalism, this being the first contribu- 
tion of this paper. With the basic ingredients at hand, 
wc present in Sec. IIIII our main result: a semiclassical 
expression for the entanglement dynamics. The formula 
is analyzed in Sec. IIVI as follows. First, the canonical 
result |16| is shown to be exactly reproduced in an ap- 
propriate limit. Second, we test the accuracy of our semi- 
classical result in describing the entanglement dynamics 
for the problem of phase coupling between two spins j. 
In Sec. |Vl wc present our final remarks. 



II. SEMICLASSICAL PROPAGATOR IN THE 
SPIN-COHERENT-STATE REPRESENTATION 

The development of semiclassical approximations for 
the quantum propagator in the coherent-state represen- 
tation has a long history. It started about 30 years 
ago with Klauder's approach [l3 on the one-dimensional 
canonzcaZ-coherent-state propagator, K(z^,z^,T) = 



-iHT/h\ 



Subsequently, other works [18|-[20| sub 



{z^\e^ 

stantially contributed to the understanding of the semi- 
classical version of K(z^,z^,r). In particular, Ref. [20| 
consists in a very detailed study of the subject and 
will be, therefore, the main support to our approach. 
Moreover, extensions of the semiclassical formula to 
further canonical degrees of freedom can be found in 
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Refs. [2l|, |22|, while derivations for spin variables are 
given in Refs. [23l - [27j . Also, it is worth mentioning a 
result on the two-dimensional semiclassical propagator 
for the case where one variable is a spin and the other is 
canonical [l^ and a recent derivation for SU(n)-coherent- 
states pgj . 

Despite this vast literature on semiclassical propaga- 
tors, only recently a result has been reported [l6|, pro- 
viding a semiclassical approximation for time-reversal 
propagators using the canonical states. In what follows, 
we extend this result by deriving a unified formula for 
the two-dimensional semiclassical propagator in the spin- 
coherent- state representation, expression which is consid- 
ered as the first contribution of this paper. However, be- 
fore presenting it, for the sake of completeness, we briefly 
review some elements of the spin-coherent-state formal- 
ism. 



A. Spin Coherent States 

Spin coherent states were introduced by Radcliffe [sO] 
in direct analogy to canonical coherent states. Since then, 
they have become important tools in a variety of areas 
of physics (see Refs. [3ll - [33 | for examples and further 
details). 

The spin coherent state associated with a particle of 
spin j is defined as 



cxp 



{sJ,} 



(1 



i2y 



J), 



(1) 



where the label s is a complex number, J-|_ is the raising 
spin operator and | — j) is the lowest eigenstate of J3 
with eigenvalue —j. The notation adopted here is such 
that both s and are dimensionless quantities. That 
is, in this paper the operator J denotes the usual angular 
momentum operator divided by h, so that its components 
satisfy 



(2) 



plus cyclic commutation relations. In terms of these 
states, an over-complete unity resolution can be written 
as 



\s){s\ dv{s) = Is, dv{s) 



2j + 1 ds(«)ds(^) 



(3) 



where s^^-* and s*-^^ are, respectively, the real and the 
imaginary parts of s, and the integral runs from —00 to 
-|-oo. In addition, spin coherent states are, in general, 
non-orthogonal as can be seen in the overlap 



|2J 



(1 + |.,P).(1+|,^J2) 



2V 



(4) 



It can be shown that \s) saturates the uncertainty rela- 
tion (Ai2)(AB2) > i|({Ai,AS})|2 [H 



for angular momentum operators, which implies that spin 
coherent states are minimum uncertainty states. 



B. Spin Semiclassical Propagator 

Let the forward — and backward = —1) quan- 
tum propagator in the spin-coherent-state representation 
be written as 



(s;,s^,T) 



(s s \e-^^"'^/^\s s ) 



Considering the limits j ^ 00 and ft — >■ with the prod- 
uct hj finite, we follow Refs. [l^, [l^, [2^ to obtain the 
semiclassical formula 



e 



(5) 



Sy) refer to different 



The indices x and y in |s) = \sx 
subsystems. We assume, for simplicity, that the spins 
have the same magnitude j, i.e., both Hilbert spaces have 
dimension 2j The right-hand side of Eq. ([S]) depends 
only on complex trajectories governed by a Hamiltonian 
function H (see below) . In terms of auxiliary variables u 
and V, the Hamilton equations are 



dH 



-2ijh Vk 



duk {l + UkVkf 



and 



dH 



2ijh ilk 



dvk {i+UkVkf 



(6) 



where k = x,y and H{u,w) ~ H{s,s*) = (s|i/|s). This 
equality implicitly defines the new variables through the 
replacement of s and s* by u and v, respectively. Trajec- 
tories contributing to Eq. ([5]) must satisfy the boundary 
conditions 



u 

u" 



and v" 
and v' 



for e 
for C 



(7) 



In our notation, single (double) prime stands for initial 
(final) time. The sum in Eq. ([5]) runs over all trajec- 
tories governed by Eqs. (j6|) and submitted to boundary 
conditions ([7|). 

The complex action = 5^(s*,s^,T) and the func- 
tion = Q^{s*,Sf^,T), in Eq. dU, are exphcitly written 



Jo 



k—x,y 



E 

k—x,y 



1 + UkVk 

diik _ dvk_ 
duk dvk 



dt + A, 



dt. 



(8) 

The factors A (accounting for the normalization) and A, 
appearing in Eqs. ([S]) and ([5]), respectively, are given by 



A^jY. ln[(l + |.^,|2)(i + |,^,|2)] 

k—x^y 

A = jJ2 ln[(l + 4«D(l+«)]- 

k—x,y 



(9) 
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At last, the prefactor of Eq. ([5]) can be written as 



n 

k=x,y 



2j 



where 



dsuvds* 



(10) 



(11) 



We point out that the phase of plays a role similar to 
that of the Maslov phase in the coordinate propagator. 
Because of the square root in Eq. ([5]), we must track it 
over time and add, after each complete turn, a phase — tt 
to the propagator. 

For future use, we differentiate to get 



dSf 



2jhds 



and 



1 + 



1 + u'^v'i^ 
v'k 

1 + « 



for ^ 

for ^ 

for ^ 

for C 



-1, 
-1, 

-1, 
-1. 



(12) 



(13) 



In addition, dS^/dT = -^ff(u',v') = -^if(u",v"). As 
shown in Appendix \^ Eqs. ((T^ and (|13l) allow one to 
write second derivatives of 5^ in terms of the elements of 
the stability matrix M, which is defined by 



5u" 
Sv" 



M 



Su' 



Mvu Mvv 



Su' 
(5v' 



(14) 



It follows that in terms of M the prefactor reduces to 




-1 



-1 



(15) 

which is clearly more appropriate for numerical purposes. 

At this point, it is worth mentioning why trajecto- 
ries contributing to Eq. ([5]) are complex in general. As 
pointed out after Eq. ([5]), s and s* were just replaced by 
the new variables u and v, respectively. However, if one 
simply considers that u = v*, a seemly natural assump- 
tion, one cannot generally find contributing trajectories 
to Eq. ([S])- In fact, for both values of ^, this would im- 
pose an excessive number of boundary conditions, since 
the evolution time T, and the initial (u',v') and final 
(u",v") phase space points would be completely deter- 
mined. This over-constrained problem can be circum- 
vented by introducing the aforementioned complex tra- 
jectories, which arc obtained by extending the real and 
imaginary parts of s to the complex plane. This pro- 
cedure is equivalent to assume that s* is no longer the 



complex conjugate of s. Such a maneuver, whose formal 
support is given in Ref. [l^l, justifies why s and s* are 
renamed into u and v. 

Finally, it is also important to note that if, in a given 
instant of time t, a trajectory has only non-null real co- 
ordinates, i.e., u(r) = [v(r)]*, then it will be always real. 
This can be seen as follows. If H is Hermitian, then 
(s|i7|s) = ((s|i/|s))*, implying that H can be written 
as a power series of the real and imaginary parts of s, 
with real coefficients. Rewriting Eq. (j6|) in terms of s'-'^^ 
and s^^\ one may verify that real points in phase-space, 
namely, those for which Imjs'^^''} = Imjs'^-^)} = 0, are 
allowed to possess only real phase-space velocities. It 
follows that the motion is constrained to the real phase- 
space. 



III. SEMICLASSICAL ENTANGLEMENT IN 
PURE BIPARTITE SPIN SYSTEMS 

The entanglement dynamics of a pure bipartite system 
composed of subsystems x and y can be quantified by the 
linear entropy of the reduced density matrix. 



S\in{Px) = 1 - PiPx), 



(16) 



where = Tr^^p, p = \^iT)){yj{T)\, and |^(r)) is the 
state of the system in a given instant of time T. The 
purity of the reduced density matrix px is given by 



Pip,) ^ Tv,{pl} = Tr, [[Tvyp{T)f} 



(17) 



a positive quantity lying on the interval [0, 1]. For pure 
bipartite systems P is symmetric, i.e., P{px) — P{Py)^ 
and keeps equal to unity for non-interacting systems. 

In what follows, we restrict our approach to situa- 
tions in which the initial state |V'(0)) is a product spin- 
coherent-state |so) = \sqx)® |sOy)j so that 5'lin(Pa;,y(0)) = 

0. By doing so, the matrix elements of the density oper- 
ator in the spin-coherent-state representation. 



(s,|p(T)|s, 



for a generic time-independent Hamiltonian become 
kernels in Eq. (|17p . In terms of the notation of the pre- 
vious section, these elements can be semiclassically ap- 
proached by 

(sj/5(r)|s^),e^i = X+(s;,So,T) A-_(sS,s^,T). (18) 

Plugging this expression into Eq. ([TT]) and taking the 
traces in the spin-coherent-state representation, we read- 
ily obtain a semiclassical version of the purity, 



Psc(T) 



X+(«,<),So,T) K^{sl{sx,Sy),T) 



X A'+ («,!«;), So, T) K^{sl,{wx,Wy),T) 
X dh'(sy) dviwy) du{sx) di'{wx)- 



(19) 
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As seen by Eq. ([3]), this integral spans the whole 8- 
dimensional real space composed of the real and imag- 
inary parts of Sx, Sy, Wx, and Wy. 

Now, let us consider, for simplicity, that only one tra- 
jectory contributes to each propagator. Then, the in- 
tegrand depends on four complex trajectories, each one 
contributing to its respective propagator and obeying dis- 
tinct, though correlated, boundary conditions, namely. 



u' = So and v" 



V = s5 and u" = {sx,Sy), 



u' = So and v" = (s*,w*), 
v' = Sq and u" = {wx,Wy), 



for ir+(«,s;),So,T), 
for K_{sl,{sx,Sy),T), 
for K+{{sl,wl),So,T), 
for K_{sl,{wx,Wy),T). 

(20) 

Although integral ([TOI is rather unlikely to be analyt- 
ically solved for general Hamiltonians, its structure is 
proper for the application of the saddle poi nt approxima- 
tion [11] . As carefully discussed in Ref. [l^l , it is possible 
to analytically extend integral ([T9| to a line integral over 
an 8-dimcnsional complex space, which is obtained by the 
complex extension of the real and imaginary parts of s^^ 
Sy, Wx, and Wy. This procedure is equivalent to working 
with the set (s^:, s* , Sy, s*, w^:, w*, Wy, w*) of eight inde- 
pendent complex variables. Obviously, such a change of 
variables implies that [see Eq. ([3])] 



dh'(s) 



2i -I- 1 ds(^)ds(^) 2j + l dsds* 



2ni (l + ss*) 



(21) 



In this new scenario, the first step of the saddle point 
method can be directly performed. It consists in looking 
for critical points {sx, s* , Sy, s*, uj^;, u)*, Wy, iZ)*) of the in- 
tegration variables. Neglectin g d erivatives of the terms 
and P^, as justified in Ref. [201, saddle points arc 
obtained from 



d 

dSy 

ds* 

dSx 

a 

dw* 



dwl 

d 

dwx 



L,^ + j-S+iiiu*x,rsl),so,T) 

L,^ + fs+iirx,w;),so,T) 

^^y + t'^-(^0' {'^x,Wy),T) 

+ f^S+iiw*x, s*y), So, T) 

Lw, + fS-{sl, {Wx,Wy),T) 



(22) 



= 0, 



where La^ = —2j In {1 + ak a^) , with ak assuming Sx, 
Sy, Wx or Wy. Using Eqs. (fT2|) and (fTS]) one shows that 
Eqs. imply that the four critical trajectories con- 
tributing to Eq. ([T^ should obey the following additional 



boundary conditions: 



Uy = Wt 



V 



w. 



and u" 

and v'x 

and u'x 

and v'' 



- Wx 

■ '^x ■ 
— * 
W^ 



for K+iiw*x,s;),So,T), 
for K_{s^,{sx,Sy),T), 
for K+i{s*x,w;),So,T), 
for K^{s*,{wx,Wy),T). 

(23) 

It follows from Eqs. and ((231) that the final bound- 
ary conditions of the four critical trajectories must be 
real, namely, u" = (v")*. Since this implies that these 
trajectories have to be real for every instant of time, we 
conclude that the critical set is necessarily composed of 
four real trajectories. Because of this constraint, the ini- 
tial boundary conditions of each trajectory become com- 
pletely determined and, moreover, turn out to be the 
same. Therefore, there is no other option but to con- 
sider that all critical trajectories actually correspond to 
the same solution departing from u' = Sq and v' = Sp. 
Clearly, this trajectory simultaneously satisfies Eqs. ([20)1 
and (p3|) . 

Now, expanding Eq. (|19|) up to second order around 
the four critical trajectories produces 



T 



dot Muu dct Mv 



e^^^AS^dviz), (24) 



where diy{z) — dv{sy) dv{wy) di'{sx) dv{wx) and 



(5z-^ = ( 5wx Swy Ss 



Swy 6sx Ssy Swl Ss* 



) , (25) 



with SPk = Pk — Pk- Here, /? assumes w or s, or still their 
complex conjugates, while k assumes x or y. In addition. 



T 



n 



(1 + sofc sSfe) 



(26) 



and the matrix A contains second derivatives of S± and 
Lqj.. The Gaussian integral in Eq. (|24p can be exactly 
solved, as shown in Appendix |B1 Using the result (jBl[) 
one may rewrite Eq. as 



Psc = 



r 



^{d-d')^-d"^ 



(27) 



where 



d 
d' 
d" 



det Muu dct Mvv _ 
det A det B + det C_det D ,_ 
det A' det B' + det C' det D', 



dct Muv det Mvu, 



with the auxiliary matrices 



A 
C 



A' D' 
C' B' 



10 
1 
10 
10 



10 
10 
10 
1 



M, 



M. 



(28) 
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In writing Pgc in terms of these auxiliary matrices, we 
have used the relations 



MuvMvv 



MvuM-i 



1 



det Mvv 
1 

det M„„ 



dctp -dctp' 
det B' det B 

det C_ dot A' 
-detC det A 



(29) 



which can be directly verified. Equation (|27p can be fur- 
ther simplified by noting that the determinant of matrix 
M can be written as 



so that 

p — 



detM = d - d' - d'\ 



r 



(30) 



det M det M + 2d" 



2<f_ 



1-1/2 



.(31) 



To derive the last equation we have used the result 
det M = T, whose demonstration is left to Appendix [C] 
Equation (PT|) is the main result of this paper. It cor- 
rectly reproduces two important properties of the quan- 
tum purity for pure bipartite systems. First, through 
the analysis of the elements of M one may readily veri- 
fies that d" ~ for non-interacting systems. In this case, 
Eq. dsn reduces to Psc(r) = 1 (and 5iin(T) = 0), as 
expected. Second, Eq. (j3ip is symmetric, since it is in- 
variant under the exchange of the indices x and y. This 
can be shown by direct inspection of Eq. (|C6p and the 
elements of d", 



dot A' = 

detB' = 

detC' = 

detD' = 



2ij7i du'^duy 



(1 + uX) d^S^ 
2ijh dvpv'^ 

2ijh du''du'^ 
(1 + uXO^ d% 
2ijh dv'^dv'^ 



det Mv 



det M„ 



det Mv 



(32) 



which are obtained from the last of Eqs. 
first of Eqs. dMI, combined with Eq. P^ . 



IV. ANALYSIS 



31) and the 



approximation method, is the solution of a classical struc- 
ture defined by equations of motion ^ and Hamiltonian 
H = {s\H\s). This result is in total agreement with those 
reported in Refs. (l3 - [l6j for canonical degrees of freedom 
and, to the best of our knowledge, is the first of this na- 
ture for systems of spins. 

Also noticeable is the fact that Psc docs not depend on 
h or j separately, except through H. A direct inspection 
of equations of motion ([5]) — the building blocks of M and 
hence of Psc — reveals an explicit dependence only on the 
product hj, which keeps finite in semiclassical regime. 
As a consequence, we expect our result to remain valid 
even in the strict classical limit, as defined by fi, = 0, 
J = 00, and hj finite. Moreover, one may regard this 
as a formal proof that entanglement must survive in the 
classical limit of closed pure systems. 

A careful inspection of the semiclassical propagators 
reveals that the exclusive dependence on hj derives from 
the fact that all four contributing trajectories coalesce 
to a single solution. As a consequence, contributions 
emerging from the exponentials, which contain, sepa- 
rately, terms on h and j, cancel out identically as evi- 
denced in Eq. (|24p . While this simplifies the calculation, 
since that all functions turn out to be expanded around a 
single trajectory, the validity of our result gets restricted. 
Indeed, it seems that semiclassical approaches containing 
just one contributing trajectory do not contemplate more 
complex behaviors, as oscillations and revivals, or even 
longer evolution times. Usually, such features are well 
described in semiclassic al phy sics only when more tra- 
jectories are considered [37h39| . Then we expect that, in 
general, our derivation be valid just for short values of 
time, region where just one trajectory is able to repro- 
duce quantum results. We point out that our program 
here was just to keep the standard steps of the saddle 
point method arriving at a first formula, letting improve- 
ments on the formalism to a future work. 

Finally, although the derivation o f Eg . (E |) demands the 
limit j — 7> c)0, as discussed in Refs. j25l uuL [28| this kind 
of approximation also applies for systems with spin j ~ 
1/2. Basically, it works because second order expansions, 
essence of the approximations performed, are enough to 
describe correctly the dynamics of spin-i systems. We 
then expect that Eq. (pij) may be also applied to this 
class of problems. 



In this section, further arguments are given which help 
one to access the qualities and limitations of the semi- 
classical formula (j3ip as a quantifier of entanglement dy- 
namics. 

We start by noting that Eq. (PT|) essentially contains 
correlations among elements of the stability matrix. Re- 
markably, this means that the onset of entanglement is 
exclusively determined by the stability of a trajectory 
departing from the center of |so). This trajectory, which 
is selected by rigid boundary conditions imposed by the 



A. The canonical limit 

A further interesting test for our result concerns the 
canonical limit. According to Refs. [H, H^], canonical 
coherent states can be obtained from spin coherent states 
through a contraction process, which is implemented as 
follows. Introducing scaled quantities s = z/^y2j and 
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j+ = \/2ja'' , one takes the limit j — !• od to get 



exp {za^} 



1 



V2 



\z) being the well-known canonical coherent state. In ad- 
dition, discarding terms smaller than j^^, it immediately 
follows that 



ss^ 



ds ds* 
ds ds* 



-{zz - zz ), 



exp <; -z^z^ 



(33) 



h dzdz* 



and 



det S, 



(?) 



n 

k—x,y 



2j 



dctS 



(?) 



(34) 

With these expressions, we convert the formalism pre- 
sented in the previous section to that of the canonical 
case. In addition, we should be still able to recover the 
semiclassical purity derived in Ref. [T6j . which is given 
by 



p(can) ^ ^-1/2 jyj^^^ 



(35) 



where 

£ = £' 



(det Muu det Mw - dct A det B) 

X (det Muu det Mw - det C det D) - £ 

£' = -4 (det Muu det Mw det A' det B')^ , 

£" = (dct A') ^ dct B det D - (det A' det B')^ 

+ (detB') detAdetC. 

In order to prove the equivalence between Eqs. ((27)) 
and ([35]) . we use Eqs. (|32|) in the limit considered to 
show that 



det A' det B' + det A' det B' 



det M„„ det M^ 



and 



\Je-e' 



det M,i,i det M^ 



det A' det B' -t- det C det D' 
6!' 



= det Muu det Mw 

- (det A det B det C det D) 
(det A'2 - detAdetC) 

^ _ det Muu 

(det B'2 - det B det D) 



dct Mv 



where the last equality was obtained by using the deter- 
minant of Eq. P^ . 



det M„„ det M^ 



det A'2 - detAdetC, 



dctB 



12 



detB det D. 



(36) 



Since T — > 1 in the considered limit, simple manipu- 
lations on the above expressions complete the proof of 
equivalence. 

Another interesting byproduct of our approach 
emerges by taking the canonical limit in only one of 
the subsystems. This procedure automatically adapts 
our formalism — after minor modifications on Eqs. (j33p 
and (|34p — to describe, for instance, spin-boson systems. 



B. Case study: phase coupling 

Let us consider two particles, x and y, both with spins 
j, coupled to the time- independent classical magnetic 
field B = (0, 0, i?3). The /ree Hamiltonian may be writ- 



ten as Hn 



H, 



(x) 



H, 



(y) 



where H, 



(k) 





for 



k = x,y. Suppose that the spins interact with each other 
via the coupling 



H = Xh X 



iv) 



(37) 



where A is the coupling parameter. The Heisenberg equa- 
tion ih{dH^''^ /dt) = [H^''\Ho+H] = implies that there 
is no energy exchange between the spins. This is why 
Hamiltonian (j37p is said to describe a phase coupling. 

Since the entanglement dynamics cannot be influenced 
by local terms, hereafter we work only with the interac- 
tion Hamiltonian (j37p instead of the total Hamiltonian 
H + Hq. Also, we assume that the initial state is given 
by \ip{0)) = \sox) (8) \soy), with 



|sofc) 



1 



2j 



j + uk). (38) 



Setting = (l+|soa;P)-' (l+|soj,P)^ and applying conven- 
tional techniques of the quantum formalism it is straight- 
forward to show that 



2A f2j\ f2j\ f2j 



rixj \n'^ 

2a, 



\SOv 



v 

XTS^S 



(39) 



where Sk = rik — n'^., Uk = "-fe + n'ki ^-^d the sum is over 
nx,ny,n'^, n'y , running from to oo . This result equals 
P{py(T)) since it is clearly invariant by the exchange of 
the indices x and y. 

In order to establish contact with the semiclassical re- 
sult, we compute the short-time expression for the entan- 
glement generation. By expanding the result p9|) up to 
second order in time we obtain 



SuniT) - 



V8\SQx \ \SQy\j\T 

{i + \soxm + M) 



(40) 
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As anticipated by the discussion of Sec. IIV Al we ex- 
pect this result to reproduce the canonical one under 
the parametrization sq/c = ^ofc / •\/2j followed by the limit 
J — > c». Evaluating the above expression in these terms, 
we obtain that 



lim 5iin(T) = 2|zo.|'ko«|'(Ar)2, 



(41) 



which indeed yields a result equivalent in structure to 
that obtained in Ref. [iBl for a system of two oscillators. 

To apply the semiclassical formalism to this system, we 
first find the classical Hamiltonian associated to Eq. (|57| : 



i?(u,v) = (v|Aft 



1 - UxVx 



, 1 + UxVx 

Equations of motion ([5]) result in 



1 - UyVy 

1 + UyVy 



(42) 



/ Ux\ 




/ Xx 











Uy 







A, 








Vx 









^Xx 





\Vy J 




\l 








-Xy 





/ Ux\ 


1 


Uy 




y ::: I 







(43) 



It is 



where A^ = iXj ( ifHjdk) and A„ = iXj ( 

clear that both UxVx and UyVy arc constants of motion. 
Then trajectories are readily obtained in terms of their 
initial conditions, 



Vxit)=v'^e-^^\ Vy{t) 



Uy{t) ^ u'ye^-^ 



-A„t 



(44) 



From them, and remembering that A^. ~ Xx{u'y,v'y) and 
Ay = Ay(u^,w^), the stability matrix is straightforwardly 
written as M = M1M2, where 



Ml 



2t 



( Xxe^^' 








' \ 





Xye^y* 














A,e-^-* 





V 








Xye-^y^ ) 



and 



M, = 



1 













1 






2Xyt 


1 — u'^^-u^^ 







1 


l-U'y^V'y^ 
















r"5 — TT 



1 



Then, as for this system T amounts to 1. and 



det A' 
detC 



l-uL'^v''-' ' 



detB' 
detD' 



~2u'^u'y\yt 

1 ^T.' ^ 5 



(45) 



we finally find that 

Psc(T) = [l + 2d"]-'/' 
1 - 



-I -1/2 



1 



(i-<,X')(i-";XO 

V8|soa:| \sQy\jXT 

(l + ko.P)(l + ko,P) 



(46) 



which agrees with the quantum result (|40p. 

This case study highlights the major difficulty of our 
approach: the semiclassical formula applies accurately 
only in the short-time regime. Nevertheless, this is not 
really surprising. As pointed out above, it is well-known 
that quantum phenomena can be well described semi- 
classically only via many contributing trajectories. As 
we have seen, this is not the case here. Actually, this 
turns out to be one of the next challenging question in 
the context drawn so far: How to improve the semiclassi- 
cal formula so as to correctly describe the entanglement 
dynamics for longer times? 



V. FINAL REMARKS 

In summary, this paper is concerned with autonomous 
systems of two spins j prepared in a product of spin- 
coherent-states. We looked at the entanglement dynam- 
ics as quantified by the linear entropy — or its kernel, the 
quantum purity — as a function of time. A semiclassical 
approximation for the purity was derived by replacing 
exact propagators by their semiclassical versions. The 
calculation, which employed the saddle point method to 
analytically solve the integrals, produced the semiclassi- 
cal expression pip , the main result of this paper. This 
formula allows one to express the onset of entanglement 
in terms of a classical structure, defined by a Hamiltonian 
function, equations of motions, and a set of boundary 
conditions involving the initial conditions. The semiclas- 
sical time-reversal spin-coherent-state propagator ([S]) is 
another original derivation of this work. 

The adequacy of our results was illustrated by some 
important analytical tests. First, the semiclassical purity 
was shown to be symmetric. This property, which is not 
trivially reproduced by classical entropic measures d, Q , 
indicates that our formula does capture the quantum 
essence of entanglement. Interestingly, however, the re- 
sulting structure is shown not to importantly depend on 
h or j separately. This constitutes a symptom of the fact 
that the semiclassical result should be accurate only in 
the short-time regime. Second, it was shown that the 
semiclassical purity correctly recovers the canonical re- 
sult [1^ in the large-spin limit. We concluded the tests 
with a case study which confirmed the accuracy of our 
semiclassical result in the regime of short times. 

Finally, it is worth noting that our results and conclu- 
sions are in consonance with many others reported for 
canonical degrees of freedom 0, 0, Hi [3l : especially in 
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what regards the link between entanglement dynamies 
and stability of underlying classical structures. A natu- 
ral continuation of this paper includes the improvement 
of the semiclassical formula so as to reproduce the ex- 
act entanglement dynamics in regimes of longer values of 
time. Work on this topic is now in progress. 
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Appendix A: Elements of the stability matrix and 
second derivatives of the action 



In this appendix we derive relations between elements 
of the stability matrix M, defined by Eq. (|14p . and sec- 
ond derivatives of the complex action 5^, defined by 
Eq. ([8]). We start by performing variations on both sides 
of Eqs. (in]) and Dealing first with ^ = +1, we get 



where 



A+ 
C, 



5u' 



B+ \ / Su" 
B+ ) \ 6v' 

(Al) 



c(+) 

''u'u' '^u'v" 

c(+) 



c(+) - 
^ab — 



daydbcc daydby 



with a and b assuming u' or v", and 



A, =-5 




1 


//2 















Rearranging Eq. (|A1[) . so as to write the final displace- 
ments 5u" and Sv" as a function of the initial ones Su' 
and Sv' , and comparing it with Eq. ((Ti)) lead to 





= Dh 


-'{' 


Muv 


= Dh 




Mvu 








Mvv 




S 


(+) " 

u'v" 



c(+) r 



q( + ) 



^u'v" 

B+, 



(A2) 



where A+ = s|;t"^, + A+ and C+ = sltl„ + 
ing these relations, one shows that 



:(+) 



''u'u' 
»u'v" 
^v"u' 

c(+) 

»v"v" 



B+M-,iMvu 



= B+M-,1, 

= D+ [Muu - MuvM-^iMvu] , 

= C+. 



Analogous relations can be found for ^ 
cntiating Eqs. (|T2]) and (fT3|) . we find that 



where 



A_ 
C_ 



Su" 

Sv' 



B_ 
D_ 



Invert- 



(A3) 



f. Differ- 



Su' 
Sv" 
(A4) 



'ab 



*u"u" 


*u"v' 


^v'u" 


»v'v' 






da^dbx 




daydby 


daydb 



with a and b now assuming u" or v', and 
A_ = -2ii, 



B 



D 



-2ijh 



f/2 














1 










1 
























1 










1 





Manipulating Eq. (jA4p in a convenient way, we get 







Sv'u" 




1 

D , 
















1 -1 ~ 












v'l 




C 


- 5 
















-1 




Mvu 


= B: 


:iA_ 




»v'u" 




Mvv 


= Bl 




"I' - J 


\_ 


v'u" 



(A5) 



where A_ 
ing them leads to 



Sy„„+A_ andC_ 



S^,^, +C_. Invert- 



s' 

*u "u ' 
»u"v' 



B MvuM-1- A , 
B_ [Mvv - MvuM„^Muv] 



D M„^, 



S[-1 = -D_M-iM„v-C_ 



(A6) 



Equations (|A2j) . (|A3)) . (|A5|) . and (|A6)) establish the in- 
tended connection between elements of the stability ma- 
trix and second derivatives of the action. In particular. 
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they prove the equivalence between Eqs. pUj) and (IT5]) . where 
provided that we identify det Ss~''g. and det g. with 



det Sitl, and det si~l 



detS^,,'^,, respectively. 



Appendix B: Gaussian integral 

In this appendix we solve the Gaussian integral 



ai = 



02 = 



14 


' det Mvu det Muu det Muv 


det M. 


hi 


Mvu Muu h:. hi Muv Mvv 




h^ 


Mvu Muu h„ Muv Mvv 


hy, 


hi 


Mvu Muu ^ Muv Mvv 


hx 




Mvu Muu ho. hi Muv Mvv 


hy, 



dl'{z), 



which the semiclassical purity Pgc depends on, as shown 
in Eq. ([M)) . While 5z and di'{z) are defined in the main 
text, the 8x8 matrix A is composed of the following 
4x4 blocks 



All 



B,7, , 



B, 



^ ' h v"v" 



with hi = (1,0) and h^ = (0,1). Using Eq. ^ and 
^■^ji « 1, which becomes exact in the limit considered, 
we finally find that 



(2j + l)4(detA) 



-1/2 



(1 



(Bl) 



B,-„, 
B... 



Bg 

( ) 



and A21 



where Cn 



-2j 



Ba, = ( (i+a.aj)^ 



ar.Ba , 



and Ba 



Bg^ 
B,,„ 





-2j 



(i+a»a;)" 



with a assuming s, s*, w, and w*. In Appendix El second 
derivatives of the actions S± are written in terms of the 
stability matrix M of the pertinent trajectory. Using 
Eqs. (jA3p and (|A6|) . and recalling that the trajectory 
associated to S+ is identical to that associated to S- , we 
rewrite the above matrices as 



A - A - ( -SMvuMu^ 

A11-A22-I 



S 

-SMuvMvv '' 



A12 = A21 
where 



SR^ 
SR, 



1 




-2i 




R,, 




1 



and 





-2i 



Appendix C: Determinant of the stability matrix 

Here we derive an expression for the determinant of M 
[Eq. (|14p ]. the stability matrix associated to the classical 
trajectory involved in the calculation of Psc- Because of 
the symplectic structure of canonical Hamilton's Equa- 
tions, the determinant of the stability matrix is constant 
and equals to 1 (see, for instance, Ref. [1^). However, 
for the spin equations of motion ([5]) , the above no longer 
holds. Our strategy to compute det M consists in in- 
troducing a new set of canonical variables Qx, Px, % and 



Py 



ml, \M, for which det Mc 



1. Then, from the 



relation between the two set of variables, det M can be 
determined. 

Assuming that Uk = Ukiqk,Pk) and Vk = Vk{qk,Pk), for 
k ~ x,y, implies that 

5w = T Sr, (CI) 

where we have defined Sw'^ = (Sux 5uy 5vx Svy) and 
Sr'^ = {Sqx 6qy 6px Spy). Non-null elements of T are 
given by the relations 



22 



OUy 

dqy 



T 9Py_ 



tl3 
t24 



hi 

ti2 



dpy 



J X 



J 9qy_ 



OUx _ _ 7 <JPx J- — 0-»x _ 7 



(C2) 



ggx 

dVy_ 

dqy 



J dpy_ 

'■^VdUy ' 



^44 



dVy 



J, 



C?ttx 
dqy 



With these arrangements, the determinant of A can be 
straightforwardly calculated, resulting that 



det A = (det S)^ K ^ ^2] 



where Jx = ^11^33 - ^13*31, Jy = ^22*44 - ^24*42, and 
the last term of each equation is obtained by inverting 
Eq. (ICTI) . 

By demanding qk and pk to be canonical coordinates. 
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one must require that 

dqk . dqk . (1 + UfcWfc)^ 



Qk = -^Uk + -T^Vk „. .. 

ouk ovk Ztjn 

{1+UkVk)^ 



{qk,H} 

an 9H 



{qk,Pk}u,..v, - g^^, 



2ijh 

dpk . , dpk . (1 + UkVk)'^ 

Pk = ^ — Uk + -^Vk = — — 

ouk ovk 2.1 J n 

{l+UkVkf 



{Pk,H}u^^. 



2ijh ^'^"'^"^'"^'"'Wk^'Wk' 

(C3) 

where Eq. ([5]) was used to ehminate the time derivative. 
In these relations, HiqxTqy^PxyPy) amounts to 



Last equahties of Eqs. (|C3p imply that 

{qk,Pk}uk,vk = Jk^ = 2^^/(1 + UkVk)'^. (C4) 

Since the stability matrix Mean in the new set of variables 
is defined by 



Sr" = Meaner', 



(C5) 



one can use Eq. (lCT|) to find that M = T" Mean (T') \ 
It follows that 

Aci T" T" I" 

det M = ^Hi4_ det Mean = ^ = r, (C6) 



H[Ux{qx,Px),Uy{qy,Py),Vx{qx,Px),Vy{qy,Py)]. 



where T is given by Eq. 
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